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We calculate Landau gauge ghost and gluon propagators in pure SU(2) lattice gauge theory 
in two, three and four dimensions. The gauge fixing method we use, sc. stochastic quantisation, 
serves as a viable alternative approach to standard gauge fixing algorithms. We also investigate the 
spectrum of the Faddeev-Popov operator. At insufficiently accurate gauge fixing, we find evidence 
that stochastic quantisation samples configurations close to the Gribov horizon. Standard gauge 
fixing does so only at specific parameters; otherwise, there is a clear difference. However, this 
difference disappears if the gauge is fixed to sufficient accuracy. In this case, we confirm previous 
lattice results for the gluon and ghost propagator in two, three and four dimensions. 
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I. INTRODUCTION 

Confinement is one of the key properties of QCD, but 
its explanation from first principles still remains to be 
found. While a linearly rising heavy quark- antiquark po- 
tential is readily simulated on the lattice, evidence for 
or against the candidates for the underlying confinement 
mechanism is more elusive. Some prominent confine- 
ment scenarios involve topological defects, like magnetic 
monopoles or center vortices, and are best described in 
specific gauges amiable towards the definition of the re- 
lated defects. It is this gauge dependence which makes 
it so difficult to pin down the confinement mechanism. 
There is some evidence that such scenarios are compat- 
ible with approaches like the Gribov-Zwanziger scenario 
IE S 9 > according to which configurations in the vicin- 
ity of the Gribov horizon account for confinement. This 
scenario makes specific predictions for QCD's Green's 
functions, e.g. in Landau gauge that the gluon propa- 
gator vanishes in the infrared (IR), the horizon condi- 
tion. In turn the ghost dressing function is IR diver- 
gent, called ghost enhancement. The Kugo-Ojima sce- 
nario [i| has similar implications. It aims at describing 
confinement by the BRST quartet mechanism and re- 
quires well-defined global BRST charges. Both scenarios 
put restrictions on the infrared behaviour of ghost and 
gluon propagators if confinement is present in the theory. 
In turn it has been shown recently that ghost and gluon 
propagators with a sufficiently large IR suppression of the 
gluon and an IR enhancement of the ghost lead to quark 
confinement indicated by center symmetry [f|. The re- 
lated IR suppression of the gluon propagator has been 
seen in both, lattice studies and continuum approaches, 
see e.g. @,0|. 

The quantities relevant for such confinement scenarios, 
like the aforementioned Green's functions, depend on the 
chosen gauge. Hence these investigations are complicated 
by the Gribov problem. This problem relates to the fact 
that local gauge conditions do not single out a unique 



gauge copy [l|, [1] , but allow for an abundance of Gribov 
copies satisfying the gauge condition on the same gauge 
orbit. 

Interestingly, Landau gauge seems not to fix uniquely 
the infrared behaviour of the propagators. In continuum 
computations with functional methods such as Dyson- 
Schwinger equations (DSE) and functional RG flows one 
observes a one-parameter family of solutions. The free 
parameter can be linked to the zero momentum value 
of the gluon propagator in the infinite volume limit, see 
[f|. It is suggestive to relate this parameter to global 
properties of the gauge fixing as Landau gauge still has 
remnant gauge degrees of freedoms. In the continuum, 
this infrared boundary condition is implemented by spec- 
ifying a non-perturbative renormalisation condition, see 
[6| . It has been demonstrated recently that a correspond- 
ing freedom exists also on a finite lattice, and that it is 
directly connected to the non-perturbative completion of 
the Landau gauge condition, [9(. Indications for this free- 
dom have already been seen in the strong coupling limit 
[3 EE, E2] where the effect is rather dramatic 13]. This 
high sensitivity of the results to the global details of the 
gauge fixing, in particular for the ghost propagator, calls 
for a systematic evaluation of the details of the thermo- 
dynamical limit of Landau gauge fixings on the lattice. In 
particular the question of whether the Gribov-Zwanziger 
scenario is indeed implemented within specific gauge fix- 
ings is of importance. 

In the present work, we compute ghost and gluon prop- 
agators in Landau gauge Yang-Mills theory on the lat- 
tice within stochastic quantisation [3, EH, EH, an al- 
ternative to standard Monte Carlo simulations. It has 
recently also received renewed interest e.g. for simu- 
lations at finite density [13, EH or in real time 
Here, we employ stochastic quantisation for gauge fixing 
[H El HI [2J, M, EE EH El IH ■ In such an approach, a 
gauge fixing force drives the Langevin evolution towards 
the first Gribov region. This potentially opens new pos- 
sibilities to circumvent the Gribov problem. It also fur- 
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thers the insights and allows for a better understanding 
of the underlying structures. At insufficiently accurate 
gauge fixing, we find evidence that stochastic quantisa- 
tion samples configurations close to the Gribov horizon. 
Standard gauge fixing does so only for specific parameter 
choices; otherwise, there is a clear difference. While this 
difference disappears if the gauge is fixed to sufficient ac- 
curacy, the sampling close to the Gribov horizon makes it 
unlikely that stochastic gauge fixing systematically sup- 
presses contributions from the Gribov horizon. Such a 
suppression would have the potential of invalidating the 
argument that in the thcrmodynamical limit only con- 
figurations at the Gribov horizon matter. In summary, 
our analysis provides strong evidence against one of the 
reasons for a failure of lattice gauge fixings to reproduce 
the Gribov- Zwanzigcr scenario. 

In the second section we define ghost and gluon prop- 
agators and discuss their general scaling properties. 
In the third section we introduce our implementation 
for stochastic quantisation and stochastic gauge fixing. 
Some of the intricacies of such an approach are then il- 
lustrated within a simple finite-dimensional model in the 
fourth section. In the fifth section we present and dis- 
cuss our results on ghost and gluon propagators in two, 
three and four dimensions, which we supplement in the 
subsequent section with some results on the spectrum 
of the Faddeev-Popov operator. We close with a short 
summary and discussion of our results. 

II. INFRARED BEHAVIOUR OF THE GHOST 
AND GLUON PROPAGATOR 



family, and the corresponding propagators obey power 
laws in the infrared, 

Urn D gl/gh (g 2 ) ~ (4) 

In d = 4 dimensions, the exponents satisfy the scaling 
relation 

k a = —2kc' (5) 

It stems from non-renormalisation of the ghost-gluon- 
vertex (65j . which has been confirmed on the lattice 
[Hi, Ei} as well as from DSE [H]. When varying d, © 
generalises to 

k a = -2k c + (d- 4)/2. (6) 

Assuming a bare ghost-gluon vertex, the predicted values 
for Kc-, which is also referred to simply as k, are k = 
i (93 - V1201) as 0.6 in d = 4, k as 0.4 in d = 3 and 
k = 0.2 in d = 2 [p]. 

The aforementioned decoupling solutions imply that 
both the gluon propagator and the ghost dressing func- 
tion are finite in the infrared. In contrast to the scaling 
solution, the decoupling solutions seem to be incompat- 
ible with global BRST invariance @, see @, EI El for 
progress towards a lattice BRST formulation. Neverthe- 
less, both solutions exhibit confinement ,as both satisfy 
the confinement criterion put forward in [5J] and both lead 
to positivity violation for the gluon. Positivity violation 
has been explicitly confirmed on the lattice in d = 3 [7l| 
and d = 4 [zl,El and also from DSE 



We are interested in the infrared behaviour of Green's 
functions in Landau gauge Yang-Mills theory. Landau 
gauge is defined via the condition, 

MS = °- (!) 

which is subject to the Gribov ambiguity. In Landau 
gauge, the gluon and ghost propagator have the form 

(D*)« = 6"» (V - D ^), ( 2 ) 

(D gh r b = -6 ab D gh ( q 2), (3) 

they are parametrised by the two scalar functions 
-Dgi/gh(<? 2 )- The standard dressing functions are given 
by q 2 Up-i /p-hfc? 2 ). Functional methods in the continuum 
such as DSE MJsOJS El, S El EI El, stochastic 
quantisation 131. 1 3 71 3811 and the functional renormalisa- 
tion group [(j, l39l. |40|. l4ll. l42l. |43| allow for a scaling solu- 
tion which is unique |4 4| , | 4 a. |46j and in line with global 
BRST-symmetry |6j, 17J, |47| . As mentioned in the intro- 
duction, there are strong arguments for a one-parameter 
family of decouplin g so lutions where scaling is violated, 

see (a IISIE! Ei, Hi El, Ei Ei, EI Ei, Ei Ei, Ei, Ei, 

EH, I62l.l63l.l64l] for related work. The unique scaling solu- 
tion is indeed the scaling endpoint of this one-parameter 



A. Lattice simulations in d = 2, 3, 4 

Lattice simulations have so far found no clear evidence 
of a scaling behaviour, with the possible exception of 
the two-dimensional case [zl Ei], which lacks dynam- 
ics [77||. Very recently it has been shown that the scal- 
ing behaviour might also be absent in two dimensions 
for the lattice Landau gauge fixings used so far [l3|]. In 
four dimensions, lattice simulations of SU(2) or SU(3) 
[Zl, Ei, El, El, Ei, EH, El, E3 yield an IR finite gluon 
propagator even on volumes up to (27 fm) 4 . DSE stud- 
ies on a torus suggest that this volume could suffice for 
the correct infinite volume behaviour to be visible [84| . 
While dissenting results concerning IR finiteness exist 
[H Ei, E3, EU , they tend to yield an IR exponent closer 
to k = 0.5, implying an IR finite propagator, than to the 
predicted value of as 0.6. 

Also for d = 3, lattice results clearly favour decoupling, 
e.g. [H, even at (85 fm) 3 [H. Recently, following a dis- 
cussion of non-perturbative ambiguities of the Landau 
gauge condition (gl, evidence has been found that the 
IR behaviour especially of the ghost propagator might 
strongly depend on the choice of the Gribov copy on the 
residual gauge orbit d. This corresponds to the choice 
of boundary conditions in functional methods [(| . 
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In the strong coupling limit [lCt I 111. 1 12, 1 9 il l , si mulations 
have likewise also been performed for d < 4 [U, [H . The 
interpretation of the results is still under debate. It has 
been argued that the lattice data are compatible with a 
scaling branch for a subset of lattice momenta [l(| , but 
also that they definitely disprove such a behaviour 
Very recent additional results provide strong evidence for 
the severity of the Gribov problem, especially regarding 
the ghost propagator [l3|. 



B. Gauge fixing problem on the lattice 

Since the ghost and gluon propagator are gauge- 
dependent quantities, their calculation requires to sam- 
ple configurations from the gauge fixing surface T = 
{A /1 (x)\d ll A ll (x) = 0}. Fixing to lattice Landau gauge 
amounts to maximising a gauge fixing functional whose 
Hessian is the Faddeev- Popov operator — d fJi D'^ b [A\, 
where is the covariant derivative. Thus, after gauge 
fixing, all configurations are located inside the Gribov 
region SI = {A^(x)\ - dD[A] > 0}. However, SI still 
contains gauge copies, in contrast to the fundamental 
modular region A C fi. On the lattice, reaching A is 
equivalent to obtaining the global maximum of the gauge 
fixing functional 



R = 



1 



N c Vd 



^Retrt/^O). 



(7) 
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This poses an NP-hard optimisation problem, equivalent 
to finding the ground state of a spin glass [93, [94[ . 



Due to Zwanziger's gauge fixing term, the equilibrium 
configurations of this equation arc fixed to Landau gauge. 
It can be shown [2f| and has been demonstrated numer- 
ically [2l[ that this equilibrium is stable only inside SI. 

While stochastic gauge fixing cannot be expected to re- 
store global BRST invariance, it is a priori possible that 
it might sample the configuration space differently from 
standard gauge fixing. E.g., it might give us a uniform 
gauge fixing on the first Gribov region SI. Note in this 
context that the statistical arguments in favour of the 
Gribov-Zwanziger scenario only work if the gauge fixing 
does not single out some region inside the first Gribov 
region in the thcrmodynamical limit. A uniform distri- 
bution avoids such a possibility. 

On the other hand, stochastic quantisation is unlikely 
to provide us with a tool for a unique gauge fixing in the 
fundamental modular region A, since it does not even 
aim at the global maximisation of the gauge fixing func- 
tional R. But this fact might not harm its prospects 
for two reasons: First, even sophisticated algorithms de- 
signed to better approximate the global maximum of the 
functional R do not yield the scaling solution on the lat- 
tice HI IH, [H, S Second, the FMR and the Gri- 
bov region should be equivalent in the thermodynamical 
limit according to a conjecture by Zwanziger [3j. So, if 
stochastic quantisation leads to a uniformly covered Gri- 
bov region, this might be already sufficient. 



B. Lattice formulation 

We simulate quenched SU (2) with the group elements 
conveniently parametrised as 



III. STOCHASTIC QUANTISATION ON THE 
LATTICE 

A. Stochastic gauge fixing 

As already discussed in the introduction, there is quite 
some evidence that Landau gauge fixing on the lattice 
still severely suffers from the Gribov problem. This al- 
ready motivates the choice of an alternative lattice gauge 
fixing algorithm. 

Stochastic quantisation is a method introduced in a 
more general context in Within stochastic quanti- 
sation gauge fixing can be implemented via a gauge fixing 
force [2fJ. An infinitesimal gauge transformation in the 
Langevin equation for the gauge field is given by 



dA« 
dt 



SS 



YM 



6A% 



Dfv b 



(8) 



l d^Afi with a gauge fixing parameter a, 



Here, v = 

and 77^ is Gaussian white noise with the properties 
(rfc(x, t)vl(x',t')) = 26 ab S^6(t - t')S(x - x'). 



(9) 



U li (x)=u°Jx)+ia a u a Jx). 



(10) 



Unquenching effects are discussed in [79j, and SU(3) is 
treated in [si|, HH, [98|, |9{|. A comparison of SU{3) and 
SU(2) is specifically interesting as most continuum im- 
plementations with functional methods use approxima- 
tions which are insensitive to l/A^-cffccts. 

The lattice formulation of the Langevin equation, ini- 
tially without Zwanziger's drift term, is [l5| 



U^x) -> exp (ia a R^ x ) U^x), 



(11) 



with R composed of the Lie derivative of the standard 
Wilson plaquette action [l00( ] . providing the dynamical 
drift force F, and a stochastic term, 



sF a 

~ fix 



= ™%s[u]. 



(12) 
(13) 



In addition, we have implemented a random walk 
method, whose results agree with those obtained from 
Langevin evolution. In this case, the form of an individ- 
ual update is 



U ll (x)^ eX p(i<r a A a Jx))U ll (x), 



(14) 
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and A a ^(x) is determined, starting from zero, by accepting 
updates 

A;(x)^A;(x)± V (15) 

with probability 

p=l h±tanhQr7F;j) . (16) 

Regardless of the method used for the dynamic updates, 
Zwanziger's additional drift force is implemented as a lat- 
tice gau ge t ransformation, sc. U^x) — > f2(x)[/ /J (x)f2 T (x+ 

fx) with [H, H, E| 

^ =op (4 AV s)- (l7) 

Here, A a is the lattice version of the Landau gauge con- 
dition, 

d-l 

A%r) = £(^(*)-^(x-£)). (18) 

By means of this implementation, the gauge force moves 
are not only tangential to a gauge orbit, but strictly on 
the orbit. 

We interchange the updates (JTTJ) and the gauge trans- 
formations 'locally', i.e., such that a gauge transforma- 
tion of all affected link variables is performed after each 
single link update. This serves to minimize the distance 
to the gauge fixing surface during the updates, measured 
by A 2 , the average of A 2 (a;). However, A 2 does not de- 
crease monotonically when a is lowered (fig. [1]), which 
presumably is due to an 'overshooting' effect. 

For a high precision of gauge fixing, we usually amend 
an interchange of dynamical and gauge fixing steps as 
described abo ve by standard stochastic overrelaxation 
(STOR) [lOl|. While this goes, strictly speaking, be- 
yond the scheme of stochastic gauge fixing, it hardly af- 
fects the 'locality' of the latter. An alternative possibility 
would be to choose a step size small enough to render the 
'amendment' dispensable, cp. the decrease of A 2 with 
the step size in fig. [TJ But this would require so small 
a step size that the autocorrelation time would become 
unacceptably large for practical simulations. 

In contrast, we refer by the term 'stand ard g auge fix- 
ing' to a usual heat-bath thermalisation |l02| followed 
by STOR. This is a 'global' procedure, as the thermali- 
sation takes place far away from T, i.e., at large A 2 . The 
difference is sketched in fig. [2J 

C. Return cycles 

It has early been demonstrated [2l| that a start outside 
the Gribov horizon leads to 'return cycles' during which 
the distance from the gauge fixing surface temporarily 
increases, i.e., the system returns to the Gribov region by 
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FIG. 1: Illustration of the change of A 2 with the gauge fixing 
parameter a on a 48 2 lattice for two different random walk 
step sizes. 




FIG. 2: Qualitative illustration of a difference between 
stochastic and standard gauge fixing. The green arrows de- 
note gauge fixing steps, the black lines connecting them de- 
note dynamic updates. A similar picture appeared already in 
0. 



a path outside the gauge fixing surface. This is expected 
since only gauge-fixed configurations inside the Gribov 
region are stable equilibrium configurations. This is by 
itself not a special property of the gauge force, as it is 
shared by every algorithm that performs numerical gauge 
fixing. But a priori, the return cycles might play a role 
in bringing about a certain distribution in configuration 
space, see subscc. IVIBI 

The existence of return cycles is readily confirmed by 
the numerics. By starting the system outside the Gri- 
bov region, e.g., at sufficiently large constant gauge field, 
and measuring the lowest nontrivial eigenvalue of the 
Faddeev-Popov operator during the stochastic evolution, 
we can indeed confirm that the configurations move to- 
wards the Gribov region, in accordance with an intuitive 
picture, see fig. [3] The lowest eigenvalue Ao increases 
monotonically, while the distance from the gauge fixing 
surface reaches a maximum before decreasing again, as 
indicated by the behaviour of A 2 . 
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FIG. 3: Left (a): Qualitative illustration of a return cycle. Right (b): Numerical data. Approach to Q, during a return cycle 
with the Langevin algorithm (24 2 ). Black: Lowest nontrivial FPO eigenvalue (= distance to dfl), red: distance to gauge fixing 
surface. 



IV. STOCHASTIC QUANTISATION IN A TOY 
MODEL 



Before we present the bulk of our results on stochastic 
gauge fixing in SU (2) lattice gauge theory in sec. [V] we 
find it useful to illustrate some issues related to stochas- 
tic gauge fixing and the configuration space in QCD by 
means of a simple toy model. 



A. Motivation 

The toy model we employ is a slightly modified gen- 
eralisation of a two-dimensional model used to illustrate 
some properties of stochastic gauge fixing in [22j . Note 
that a different t oy m odel for stochastic gauge fixing has 
been studied in [l03j |. As a modification of the model 
from [22j , we have chosen the sign of the action such that 
the dynamical effects and the gauge fixing effects can be 
clearly distinguished. The generalisation consists in ex- 
tending the model to a higher number of dimensions. It 
is motivated by the possibility of observing an 'entropy 
effect' that leads to an accumulation of configurations 
close to the 'Gribov horizon' of this model, since such 
configurations should account for confinement in QCD 
according to the Gribov-Zwanzigcr scenario. Hence the 
model could give an idea about the dependence on the 
three effects: the gauge fixing force, the dynamics of the 
action and the entropy. 



B. The toy model 



where x 2 = Y^—i x 2 and the gauge fixing parameter a 
is positive. Notice that this force is not the gradient of 
an action. It drives the configurations towards the gauge 
fixing surface y = 0. The set 



n = 



y = 0,x < 1 



(20) 



may be referred to as the 'Gribov region' of this model. 
This terminology is sensible because the Gribov region is 
where all configurations are in the limit of perfect gauge 
fixing - and because the asymptotic distribution for a — > 
is 

dP(x, y) oc x n ~ l {l - x 2 )6(l - x 2 )S(y)dxdfldy, (21) 

where d£l is understood as the differential solid angle in 
n dimensions. 

The 'Faddeev-Popov determinant' of this model is 
— a _1 (l — x 2 ), which is proportional to (f2Tj) only for 
n = 1 and vanishes, like (|21[) . linearly near the horizon 
x = 1. The FP operator is obtained from the linearisa- 
tion 4iv 

at 

v = y. 



hv, with v = the gauge condition, thus here 



n-l 
n+1 ' 



The distribution (|2ip is maximal at x = 

The maximum approaches the 'Gribov horizon' dQ as 
n grows. It is interesting to sec whether stochastic gauge 
fixing supports this effect. Obviously, the gauge fixing 
force is orthogonal to the gauge fixing surface in its im- 
mediate vicinity, see fig. BJ since this surface is the equi- 
librium of the gauge force. In addition, fig. 0] illustrates 
that this equilibrium is stable only for x < 1, which jus- 
tifies the definition ([20]) of the 'Gribov region'. 



In the model, the variables Xi € K, (i € {1, . . . , n}) and 
i/eE are subject to two forces, a gauge fixing force and 
a force derived from an action. The former reads 



K 9 



f{\~x 2 )y 
2x t y 2 



(19) 



C. Numerical simulation 

To see explicitly the effects of the various parameters, 
a numerical simulation of this toy model is easily possi- 
ble, e.g. by means of a random walk algorithm, in which 
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FIG. 4: The gauge fixing force of the toy model for n 
a' 1 = 0.15. 
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FIG. 5: 'Entropy effect' for the (not gauge-invariant) distri- 
bution of values of x at the parameters a = 5 ■ 10 _J , j3 = 0, 
step size r\ = 0.001. 



the probability for a local update with fixed step size r\ 
(introducing {x„}, v = 0, . . . ,n with y = xq) is given by 

p(x v ^x v ± r)e v ) = - ± - tanh (^({a^})) . (22) 

A gauge-invariant dynamics may be defined from the or- 
bits 



= Xi exp ( — - — y ) (= const.), i = 1, . . . , n (23) 



S^(3u 2 , u 



i=l 



(24) 



5 is monotonously increasing for < x < 1 and has 
therefore the effect of driving the configurations toward 
x = 0, similarly to the expected effect of the QCD action. 
The dynamical part of the force is given by 



K a 



dS du 
du dx v ' 



(25) 
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FIG. 6: Increasing j3, the strength of the gauge action, at 
n = 500 for a = 10~ 4 , . . . , 10~ 8 . Values of f3 indicated in the 
plots. Step size r\ = 10~ 4 . 
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and the drift force of the random walk algorithm, eq. 
(HU), is 

K v = K a v + K*. (26) 

In the following figures we show the distribution of gauge- 
fixed configurations, p(x). The increasing proximity of 
the configurations to the Gribov horizon with larger n is 
illustrated in Fig. [5l The effect of the gauge fixing force 
is as expected, since the numerical results obey the distri- 
bution pTj) . The most interesting observations pertain 
the combined effect of the three parameters: n, a and 
(3. As illustrated in Fig. [51 keeping n fixed, the action 
has little effect up to some value of (3 (compare figs. [6jt 
and [6b), and then suddenly starts driving the configura- 
tions toward the interior of the Gribov zone, selectively 
depending on a, sec figs. [SJ; and [SJi. This happens at 
smaller values of (3 for lower n. 

Notice that in all cases there is a small proportion of 
configurations just outside the Gribov horizon. This is 
due to the fact that the confinement to the Gribov region 
is achieved by return cycles engendered by the unstable 
modes which K 9 develops beyond the horizon, see fig. 2] 
The fraction of configurations outside the Gribov horizon 
shrinks as the strength of the gauge fixing force increases, 
see fig. |Bk, 

The lesson from this toy model is then that there are 
subtle and non-linear effects in the combined action of 
the gauge fixing, the gauge dynamics and the entropy, 
and we should not be amazed to find such non-trivial 
behaviour in the realistic case of lattice gauge theory. 



V. RESULTS FOR THE GLUON AND GHOST 
PROPAGATOR 

A. Technicalities 

The relation between the continuum and the lattice 
momenta is, as usual, 

q ^ )= w)H^)- (27) 

We impose the cylinder cut |l04l | on the momenta. The 
propagators have been renormalised at \i = 2.5 GcV. For 
some details concerning the implementation of the ghost 
propagator, see subsec. IV CI 

We compare our results for the IR exponents with the 
'scaling predictions' [3l|, which are re = 0.2 in <f = 2, re « 
0.4 in d = 3 and k sa 0.6 in d = 4, with ka determined 
by eq. ©, seeHU 




q [GeV] 



FIG. 7: Gluon propagator on a 200 2 lattice at (3 = 10 (a = 
0.082fm). 

B. Gluon propagator 

We have calculated the scalar gluon propagator func- 
tion in two, three and four dimensions, 

with the standard lattice discretisation of the gauge field, 
A^(x) oc tt°(i). We find significant differences between 
stochastic and standard gauge fixing at relatively large 
A 2 . In turn, if we fix the gauge to a high precision, 
A 2 < 10~ 15 , the differences disappear, and we confirm 
previous lattice results in all dimensions. Notice that at 
large A 2 the picture is complicated by the fact that the 
results of standard STOR after a heat-bath thermalisa- 
tion strongly depend on the value of p, which fixes the 
probability of non-optimal gauge transformations in the 
STOR algorithm. The differences mentioned are directly 
related to the distance from d£l, which is visible from the 
spectrum of the Faddeev- Popov operator, see section IVT1 

1. Two dimensions 

The two-dimensional case is special : the re are analyti- 
cal results for, e.g., the string tension [l05j |. and the exis- 
tence of scaling behaviour has been viewed as relatively 
uncontroversial also from lattice data [75|, [7(| . Our re- 
sults for this case appear to corroborate the latter find- 
ings. The data are shown in fig. [7] We extract the 
infrared exponent ka = —1.37, which corresponds to 
K = 0.19 if eq. ([6]) is assumed to hold. This extrac- 
tion is done via a power law fit to the five most IR data 
points after discarding the two lowest non-vanishing mo- 
menta, which closely resembles the method employed in 
[75l ] . x 2 /ndf is safely below 1. As the value of k nicely 
agrees with the prediction of 0.2, this appears to speak 
in favour of the scaling solution in 2D on the lattice, cor- 
roborating previous results obtained with standard gauge 




FIG. 8: Top (a): Gluon propagator in 3D at (3 = 4.2 (a = 
0.17 fm) on 80 3 (black circles) and 150 3 (red squares). Bottom 
(b): Blowup of the same data in the IR. 



fixing |75l . 17a ]. Even more importantly, a comparison 
with the ghost propagator (subsec. IV C|) shows that the 
scaling relation ^ is approximately fulfilled. Note, how- 
ever, the surprising fact that the naive scaling hypothesis 
even fails in two dimensions at (3 = 0, see 1 1 3I J - Here it 
has been found that while approximate scaling holds true 
for the gluon propagator for intermediate momenta, the 
ghost propagator is subject to drastic changes if changing 
the global properties of the gauge fixing. 



2. Three dimensions 

In three dimensions, it has been observed in lattice 
simulations [8^ | that the gluon propagator peaks at q = 



350 



+ 100 
-50 



MeV, but still approaches a non-vanishing value 
in the infrared, D(q 2 — > 0) > 0. We confirm this obser- 
vation with stochastic gauge fixing on a 150 3 lattice at 
(3 = 4.2 («(26fm) 2 ), see fig. |S] Furthermore, we observe 
a sign of a finite volume effect when comparing the data 
to the 80 3 case (blowup i n fig . [SJa). In agreement with 
previous results for d = 3 [80|, there is no evidence that 
this suffices to yield the scaling solution as the volume 
V ^ oo. 




FIG. 9: Top (a): Gluon propagator, 40 4 , (3 = 2.2 (a = 
0.21 fm). Bottom (b): For comparison, the result of a sim- 
ulation with lower g. f. accuracy, sc. only A 2 < 10~ 8 , but 
where also the 'fine-tuning' is done by stochastic gauge fixing 
(see text). 



This is compatible with the findings in [l3j for (3 = 0. 
There it has been shown that the scaling window within 
a standard gauge fixing moves towards larger momenta. 
Naively translated to non-vanishing (3, this suggests that 
the assumed scaling region has to be washed out by the 
dynamics if we go to higher dimensions. 

For comparison, we have also performed simulations 
with standard gauge fixing on a 150 3 lattice, which yield 
a gluon propagator essentially identical with the one ob- 
tained from stochastic quantisation. 



3. Four dimensions 

For d = 4, we again confirm earlier lattice results 
[80l HH for sufficiently accurate stochastic gauge fixing. 
The gluon propagator on a 40 4 lattice at (3 = 2.2 is, like in 
3D, clearly finite in the IR. Unlike in 3D, it does not pos- 
sess a peak at finite momentum, see fig. |H1 V ~ (8.4 fm) 4 
is below the maximal lattice size already investigated 
with standard methods e.g. in [10, Hi|. But the lower- 
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dimensional results, including those on the spectrum of 
—dD, sec. I VII suggest that the problem is not a finite 
volume effect, but may rather be the well-known gauge 
fixing problem, as laid out in sections HI and HT1 

Our value of the gauge functional R = 0.8278(1) essen- 
tially agrees with the value obtained from standard gauge 
fixing, while it is significantly below the value after im- 
proved gauge fixing aimed specifically at maximising R 
[951 ]. Here, 'standard gauge fixing' denotes a method of 
choosing first copies after overrelaxation, and 'improved 
g.f.' refers to employing simulated annealing and Z2 
flips. This is consistent with our other findings, indicat- 
ing that the effect of stochastic gauge fixing is similar 
to the one of standard g.f., which renders the putative 
lattice gauge fixing problem more general. 

Since our fine-tuning by stochastic overrelaxation may 
change the location of the sampled Gribov copies, we 
have also produced some results where the fine-tuning has 
been performed entirely within the scheme of stochastic 
quantisation. In case of a random walk implementation, 
this amounts to decreasing 77 — > by some prescription, 
e.g., exponentially. Even though the accuracy of gauge 
fixing has been lessened here, the result is virtually the 
same, see fig. [2b- 



C. Ghost propagator 

The ghost propagator on the lattice is given by 

£> gh (fc) = E 



. a b 



Jk-(x-y) 



(29) 



with k ■ x = 27r^^ kfiXfi/L^. Its calculation is computa- 
tionally expensive, since it requires the inversion of the 
lattice Faddeev-Popov operator M, whose standard form 
may be found e.g. in We ha ve im plemented both 

the point source method, see e.g. [106| |. and the plane 
wave source method, see e.g. [l07j . This means that M 
is inverted either on a point source s^' x = 5 ab (6 x .o — 1 /V) 
or a vector of plane waves s%' x (k) = 5 ab e lk ' x . While the 
former method leads to heavy fluctuations in the UV, 
the latter one is much more time-consuming. The re- 
sults shown here have been produced with a plane wave 
source, but sometimes inverting the Faddeev-Popov op- 
erator only on a proper subset of the momenta surviving 
the cylinder cut. 

In two dimensions (d = 2), we extract from the ghost 
propagator, like from the gluon propagator, an infrared 
exponent close to the expected value re = 0.2, sc. n = 
0.17. The basic fitting method is the same as for the 
gluon data, see IV B ll See fig. [TOb for results on a rather 
large lattice. 

For d = 3, we expect to see a finite- volume effect for 
the ghost propagator for the volumes used here, 40 3 and 
80 3 . A fit of the 40 3 data in fig. [TO] to a power law yields 
K = 0.25, which is already clearly below ss 0.4, i.e., the 
'scaling prediction'. A fit of the 80 3 data in fig. [TO] to a 
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FIG. 10: Ghost dressing function in (a) d = 2 (top), (b) 
d — 3 (middle) and (c) d = 4 (bottom), at /3 = 2.2, 4.2 and 
10, respectively. 



power law yields k « 0.2, indicating clear finite volume 
effects. 

For d = 4, a similar effect is apparent from fig. 1 10b . 
On a small lattice, sc. 20 4 , we obtain k = 0.26. Again, 
this is far below the 'scaling prediction' of w 0.6. The 
discrepancy is significantly larger than for d = 3. More- 
over, the scaling relation ([5]) is violated thereby, given 
the gluon result of ka = —1 (subsubsec. IV B 3[) . which 
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implies k = 0.5 under ([5]). On the larger lattice, the IR 
behaviour is again harder to fit precisely, but clearly less 
divergent. Here, we have produced additional data with 
the point source method. This introduces heavy fluctu- 
ations in the UV, but does less harm in the IR, which is 
our focus of interest here. From these data, we obtain 
K = 0.19. - The fact that k doesn't v anish may well 
be due to the finite lattice extension [H, 1 108f | or also to 
choosing the first instead of the 'best' Gribov copy (97j . 



VI. SPECTRUM OF THE FADDEEV-POPOV 
OPERATOR 

A. Motivation and technicalities 

As the first Gribov region O is the set of configurations 
with positive Faddeev-Popov operator (FPO), —3D > 0, 
the lowest eigenvalue of —3D vanishes at the Gribov hori- 
zon dfl. Hence, its value may be interpreted as a measure 
of the distance from d£l. This is an interesting test of 
the Gribov-Zwanziger scenario of confinement, since this 
predicts that configurations near d£l account for confine- 
ment. The lattice discretisation of the FPO is a 3V x 3V 
real symmetric matrix. Previous lattice studies have es- 
tablished that the gauge copies tend to move towards the 
Gribov horizon as V is increased, and away fro m it as th e 
gauge fixing is improved, i.e., as R is increased [l09lllld ]. 

We want to compare the spectra after stochastic and 
standard gauge fixing. The fact that we intend to do so 
also at insufficient gauge fixing (i.e., relatively large A 2 ) 
imposes a restriction on the method used to extract the 
eigenvalues of the FPO. In particular, we cannot assume 
that the Faddeev-Popov matrix is positive definite. This 
prevents us from applying the conjugate gradient method 
for the matrix inversion and exploiting its connection to 
the Lanc zos algorithm for finding a subset of its eigenval- 
ues [llll ]. Thus, we find all eigenvalues of the 3V x 3V 
Faddeev-Popov matrix with a standard algorithm for real 
symmetric matrices, sc. symmetric bidiagonalisation and 
QR reduction, which is rather impracticable on large lat- 
tices. 

The standard lattice discretisation of the FPO which 
we use is the Hessian of the gauge fixing functional J7J). 
As usual, terms that vanish under the gauge fixing con- 
dition are omitted. To employ this discretisation also at 
imperfect gauge fixing is clearly not a unique choice, but 
has the advantage that the three zero modes of the FPO 
exist to a better approximation. The 'lowest nontrivial 
FPO eigenvalue Ao' is always understood without these 
three trivial eigenvalues. While this choice requires infer- 
ences from the value of Ao to the distance to the Gribov 
horizon to be taken with a grain of salt, it still permits a 
qualitative comparison between stochastic and standard 
gauge fixing. 



B. Results 

A few results on the Faddeev-Popov eigenvalues have 
already been presented in the analysis of return cycles, 
subsec. IIII CI the bulk of results is shown in the following. 

While we have investigated the full Faddeev-Popov op- 
erator spectrum mainly on tiny lattices, like 24 2 , the pic- 
ture we have obtained is rather unambiguous. Let us first 
state the results for stochastic gauge fixing, i.e., random 
walk steps with 'locally' included gauge fixing. At some 
intermediate A 2 , i.e., at a not very small distance from T, 
the distribution of the lowest non- vanishing eigenvalue of 
the Faddeev-Popov operator after stochastic gauge fixing 
shows a very pronounced peak close to the Gribov hori- 
zon. The peak is located slightly outside the Gribov hori- 
zon, see fig. [TT] (upper plots). This implies that Landau 
gauge is not fixed properly at this A 2 , since — dD is the 
Hessian of J d x\ a A\ 2 , which is minimised by numerical 
Landau gauge fixing. 

In the case of standard gauge fixing, i.e., STOR gauge 
fixing after heat-bath thermalisation, the distribution of 
Ao exhibits a twofold dependence on A 2 and the STOR 
parameter p. It resembles its counterpart from stochas- 
tic gauge fixing only around some specific value of p. Let 
us spell this out in more detail. At intermediate A 2 and 
large p, the procedure yields a peak at some distance 
from the horizon and clearly inside of it. When the ac- 
curacy of the standard gauge fixing is further lowered by 
stopping STOR at a larger A 2 , the distribution is shifted 
outside of the horizon, without a peak especially close 
to d£l. At intermediate A 2 , we compare different val- 
ues of p. At small p, much more configurations are out- 
side the horizon, and at intermediate p, the observation 
from stochastic gauge fixing is qualitatively reproduced, 
see fig. [TT] (lower plots). Interestingly, the latter value 
roughly minimises the number of gauge fixing sweeps to 
obtain the desired value of A 2 . For sufficiently small 
A 2 , the distribution obtained from stochastic gauge fix- 
ing coincides with the one from standard gauge fixing 
regardless of p, and no bias towards the Gribov horizon 
remains. This is evident from the contrast between the 
left plots of fig. [TTIand their counterparts at smaller A 2 , 

fig.ni 

Moreover, fig. [T2] illustrates that at the same inter- 
mediate value of A 2 , the gluon propagator is less similar 
to the fully gauge-fixed result if more configurations are 
outside the Gribov region (Ao < 0). This difference is 
created, as described, by varying p. The gluon propaga- 
tor at vanishing momentum is significantly anticorrelated 
with Ao, see fig. [13] 

A peak close to <9f2, but of course inside f2 is of inter- 
est especially for the Gribov-Zwanziger scenario of con- 
finement, since it would imply that those configurations 
are preferably sampled which account for confinement 
according to this scenario. However, generating config- 
urations at smaller A 2 by adjusting a and the step size 
does not help to generate such a peak: For small A 2 , the 
distributions from stochastic and standard gauge fixing 
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FIG. 11: Histograms of the lowest nontrivial FPO eigenvalue Ao in d = 2 resp. d — 3 from stochastic (top) and standard gauge 
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FIG. 12: Gluon propagators corresponding to the lower left 
plot of fig. [TT] 



become indistinguishable, regardless of p. This is also 
true if A 2 is lowered gradually by decreasing a and e 
rcsp. rj simultaneously. 

It is possible to relate the distribution after stochastic 
gauge fixing at intermediate A 2 to a speculative scenario 
involving 'return cycles', see subsec. IIII CI In this in- 




FIG. 13: Scatter plot of lowest FPO eigenvalue vs. gluon 
propagator at q = for intermediate A 2 (data like in lower 
left plot of fig. Ill [) . Correlation coefficient: p = —0.40 (p = 
0.1, N = 36000 configurations) resp. p = -0.29 (p = 0.97, 
N = 48000). For clarity, only 10% of data points are shown. 



tuitive picture, the Langevin evolution leads to an accu- 
mulation of the gauge-fixed configurations at the Gribov 
horizon. 



12 



600 



400 



2oo 



stochastic gauge fixing, A -1.5- 10 

A 2 <10- 10 

A 2 <10-' 5 





le-05 



0.0001 



-0.02 -0.01 



0.01 



0.02 
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2 < 1CT 15 ). Bottom (b): Standard gauge fixing 
for different values of p. 
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FIG. 15: Scatter plot of ghost propagator at lowest non- 
vanishing momentum vs. approximate lowest nontrivial FPO 
eigenvalue. Top (a): d = 2. Correlation coefficient r = —0.82 
resp. —0.83 in all three cases. Bottom (b): d = 3. r = —0.74 
resp. —0.79 in all three cases. Between 50 and 100 configura- 
tions per lattice size. 



Moreover, we find a strong anticorrelation between the 
IR ghost propagato r and the lowest FPO eigenvalue, con- 
firming findings of [l09l . Ill61 | with standard gauge fixing; 
sec fig. [T51 f unrcnormafiscd data). Note that on the larger 
lattices, we have not calculated a ll FP O eigenvalues, but 
rather employed the method of jlll| . Thus, in fig. [151 
'Ao' does not refers to the strictly lowest eigenvalue, but 
to the smallest one as found by the Lanczos procedure 
related to the conjugate gradient method. 

Finally, we also observe a finite volume effect, sc. that 
the lowest eigenvalue tends to be s maller on larger lat- 
tices, again like already reported in [l09l . Ill0| . 



VII. CONCLUSION 

In the present work we have studied lattice Yang-Mills 
theory in Landau gauge within stochastic quantisation 
in two, three and four dimensions. We also have stud- 
ied a finite-dimensional toy model in order to illustrate 



some subtleties of gauge-fixed Langevin evolutions. In 
this model wc see an interplay between the dimension- 
ality of the system, the gauge force and the dynamics 
induced by the action. 

Our results for the ghost and gluon propagator and 
the low-lying eigenvalues of the Faddeev-Popov opera- 
tor agree with the lattice results obtained so far with 
standard gauge fixing procedures. The detailed analysis 
revealed an interesting approach towards the final distri- 
bution on the first Gribov region if the gauge fixing is 
bettered. For worse gauge fixing quality the distribution 
shows also a sharp peak slightly outside the Gribov hori- 
zon. This implies a statistical bias towards the Gribov- 
Zwanziger scenario. We conclude that stochastic quanti- 
sation provides a tool of how to avoid a possible mech- 
anism against the Gribov-Zwanzigcr scenario. Still this 
does not suffice to adjust the gauge fixing appropriately, 
i.e., to obtain on the lattice the full one-parameter family 
of solutions found in continuum studies. This could sug- 
gest that the gauge fixing problem on the lattice might 
be even more general than previously assumed. 
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Together with further evidence from /3 = computa- 
tion on the lattice [l3[ , as well as recent work on Landau 
gauge fixings Q , our results hint at the severity of the 
Gribov problem in lattice simulations. We hope that this 
problem can be resolved in near future al ong the ways 
outlined in [l3[ , and in particular in [gji Ill2j . 
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